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A comparative analysis has been done of the formerly established two self-consistent solutions for 
the density of quasiparticle states in doped d-wave superconductors, having strikingly different and 
disputed behavior near the Fermi energy. One of them (1) remains finite in this limit, while the 
other (2) tends to zero. To resolve this discrepancy, the known loffe-Regel criterion for band states, 
widely used for doped semiconductors, was applied to these solutions. It is shown that both them 
are valid in limited and different energy regions, where the corresponding quasiparticles are weakly 
damped. In particular, density of states of nodal quasiparticles near the Fermi level is provided by 
the (2) solution, while the (1) only applies far enough from this level. 
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The self-consistent T-matrix approximation (SCTMA 
[11, or FLEX method Jj) is extensively used for descrip- 
tion of quasiparticles and their density of states (DOS) in 
disordered crystals, in particular, in the doped high-Tc 
superconducting cuprates. Its advantage consists in rela- 
tive simplicity and transparency for numeric calculations. 
However, this simplicity can be sometimes mischievous if 
the solutions are used without limitations on their va- 
lidity. In fact, a number of controversies emerge when 
comparing some SCTMA results with those of other 
approaches 0,0, which suggest the need for a more de- 
tailed substantiation of the method. Perhaps the central 
point in this discussion is now the question about the 
quasiparticle DOS p(e) close to the center of supercon- 
ducting gap, that is at £ ^ 0, where the DOS of clean 
d-wave superconductor vanishes as Pci(e) oc const • |e|. 
The theoretical predictions for its behavior in presence 
of impurity scattering include: i) various kinds of ten- 
dency to zero 3, 13, 0, 0, a) tendency to a constant 
value 0)0) ^'^d even in) to infinity Such extreme 
ambiguity motivates us to reconsider this situation in a 
more general context of the theory of elementary excita- 
tions in disordered systems fll| . 

There is a consensus on that the single-particle exci- 
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FIG. 1: General Mott picture of spectrum of elementary ex- 
citations in a disordered system. 



tations in disordered sytems can be either of band (ex- 
tended) type or localized type ^3|) both types forming 
certain continuous regions of spectrum separated by the 
so-called mobility edges £c (Fig. The extended states 
can be approximately described by the wavevector k, 
through a dispersion law and a broadening Ek, as 
far as the loffe-Regel criterion (IRC) fl3( for the mean 
free path i and the wavelength A holds: 



£> A, 



or fc|9£:k/ak| >Ek. 



(1) 



Then the real and imagnary parts of self-energy Sk 
in the disorder averaged Green function (GE) Gk = 
(e — Ek — Sk)~^ permit to define: -Ek = Ek — ReSk (Ek), 
Ek = ImEk(Ek) Hi, and the modified DOS p{£) = 
TT^^ImG « /Oo(e)/ (9Ek/9ek)£;^=£ (^k and po(e) being 
respectively the dispersion law and DOS in pure crystal). 
The condition £ ~ A is reached when Ek approaches the 
mobility edges, then the very notion of self-energy cor- 
rection to the initial band spectrum ek ceases to make 
sense and the averaged properties of localized states are 
only described by their DOS p(e). 

In the case of disorder due to fixed impurity perturba- 
tion Vl at random sites p (Lifshitz model) for the normal 
metal quasiparticles, the modification of band spectrum 
can involve new specific features like local or resonance 
levels 15j, and there are various ways to expand the 
self-energy in groups of interacting impurity centers \\(\ , 
analogous to the classical Ursell-Mayer group expansions 
(GE) for statistical su m llTl . Eor instance, the so-called 
fully rcnormalized GE [l8j reads: 
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where A 



p,p' 



T4Gp^pV(l - GVi) includes Gp,p' 



^"^EkVke"''^^"^'''Gk and G = Gp,p. Other types 
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of GE can differ from Eq|5] either in the structure of 
next to unity terms and in the degree of renormaUza- 
tion of G and A functions present in them. The rel- 
evant expansion parameter is not simply the impurity 
concentration c = Sp-^~^ (supposedly small, c 1) 
but the "gas parameter" cX^nT^o^o.n for the "non-ideal 
gas" of impurities with effective interaction described by 
the (energy dependent) functions Ap pi. Hence the con- 
vergence of the series (0) turns also energy dependent, 
reflecting the division between the above referred types 
of states. It can be shown that this convergence is equiv- 
alent to validity of the IRC Within the energy do- 
main of convergence, the self-energy can be approximated 
by Eq|2] with only unity term retained in the brackets: 
Ek Ti — cVl/(1 — GVi,), the momentum- independent 
SCTMA form. But beyond this domain, the SCTMA 
does not make sense and a better description of DOS is 
obtained with GE's, different from Eq[21 Below we check 
the fulfillment of IRC for nodal quasiparticle states in 
a d-wave superconductor with dopants and conclude on 
the validity of the SCTMA solutions. 

Let us start from the most common model Hamiltonian 
for this problem: 



'(''-''>^t^3V'k', (3) 



N ^ 



k,k,,p 



Here the Nambu spinors Vk ~ {A^^^^-^-ij include 
Fermi operators of normal quasiparticles with the sim- 
plest 2D dispersion law = 2t{2 — cos akx — cos aky)—fi 
in square lattice, approximately constant DOS p{e) « 
po — 4/ {ttW) where W = 8tis the bandwidth, and chem- 
ical potential /i; tj are the Pauli matrices. The d-wave 
gap function is Ak = AO (e|j — cos2(y9k, where the 
"Debye" energy ed ^ A, and (pk — a,rcta,n ky / defines 
the nodal lines k^ — ±ky. The Lifshitz perturbation term 
in Eql^l produces scattering of quasiparticles, modelling 
the impurity effect of dopants. 

The relevant GF is a Nambu matrix Gk = (^(^V'klV'k^^ 
with matrix elements being Fourier transformed two-time 
GF's: 



All the impurity effects are now accounted for by a GE 
for the self-energy matrix Ek |^ (cf. Eql^J: 



J2 [Ap,p'e^'^(p'-p)-Ap,p,V,p 



P't^P 



X 1 ^p,p'^p',p 



(5) 



Here the matrices: V — Vl%, G ^ N ^ Gk, and 



-1 ^ e»k'(p-p')gj^, (i-Gf) \ 



k'5^k 



and some additional restrictions are imposed on sum- 
mation in momenta in the products like ^p,p'^p',p, re- 
sulting from a specific procedure of consecutive elimina- 
tion of GF's in the infinite chain of coupled Dyson equa- 
tions. There are possible different such procedures and, 
respectively, different types of GE 0. Generally, GE's 
are only asymptotically convergent and the best choice 
between them is determined by their convergence range 
with respect to energy e. 

The conditions for convergence of different GE's were 
studied in detail for a number of types of elementary 
excitations in crystals with impurities [131,13; and this 
permitted to establish certain general criteria for the cor- 
responding characteristic regions of spectrum. In partic- 
ular, the region of band states is best described by the 
so-called fully renormalized GE which ceases to converge 
at approaching the mobility edges where IRC, Eq^ fails. 
For the GE, this is expressed by the tendency of all its 
terms, next to the unity in curled brackets in Eq[3 to 
become ~1. 

Alike the above mentioned scalar case, the SCTMA 
just corresponds to the fully renormalized GE, restricted 
to only its first term. Hence it is only justified when IRC 
holds. Bearing this in mind, let us analyze the SCTMA 
solutions in the vicinity of the Fermi energy for the sys- 
tem, described by the Hamiltonian 01 

Then, using the k-independent SCTMA self-energy 

Yi = V {l — GV^ and following the procedure of Refs. 
[2^,0, one can arrive at the explicit average local GF: 



{{a\h)),^i / e^(^+«)*({a(i),6(0)})dt 



where (. . .) is the quantum statistical average with the 
HamiltonianOand {., .} is the anticommutator of Heisen- 
berg operators. In analogy with the above scalar GF Gk 
for normal quasiparticles, the general solution for this 
matrix is 



Gk = - M - CkT'a - AkTi - Ej 



(4) 



G = -TrG = epo 



=K 



(6) 



In EqlBJ the renormalized energy e — e — E (e) in- 
cludes the scalar value E — TrE/2, the parameter = 
— p.po/2), and K is the 1st kind full elliptic integral. 
Having the explicit relation E = cV^G/{l - V^G^) 
between the self-energy and GF, where the renormalized 
perturbation parameter V = -I- Vi,pQ ln(l — Jl/p)], 

one obtains the self-consistent equation for G which can 
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FIG. 2: Real and imaginary parts of self-consistent Green 
functions satisfying Eq|S| Solid lines are for G*-^-* solution 
with the asymptotics, EqQ (dashed line), and dash-dotted 
lines for G^^' solution. 



be solved in principle numerically. Since the analytic 
structure of Eq|Sl involves singular points in the com- 
plex G plane (including essential singularity of the K- 
function), it possesses multiple solutions. The physical 
solutions among them should be then selected by IRC, 
as a necessary condition for SCTMA validity. 

The analysis turns most transparent in the important 
limit e — !■ (related to the Fermi energy). There are 
two characteristic solutions of Eq|Sl in this limit. 
One of them, G (e) = G^^^ (e), tending to a constant 
imaginary value, G(i) (e -> 0) i ■ const, was first ob- 
tained by Gor'kov and Kalugin in the Born scattering 
limit 's'l and then by P.A. Lee in the unitary scattering 
limit 9^. Later on, it was repeatedly reproduced by var- 
ious numerical techniques [5j and hence believed to be 
the unique SCTMA solution. Neverthless, it was shown 
recently by the authors 7] that another solution exists, 
G (e) = G*^^' (e), with low energy asymptotics: 



G(2) (e) « ^ 



A 



TTcV^po In 



, (7) 



which tends to zero with e. The behavior of real and 
imaginary parts of the two solutions in function of energy 
for a particular choice of parameters: c = 10%, Vl = 0.3 
eV, W ^ 2 eV, and A = 30 meV, is shown in FigEl 
Notice that at low energies, e A, the solution G^^^ is 
dominated by the above mentioned imaginary constant, 
presented as incpogo where 170 < 1 is a root of a certain 
transcendental equation . In contrary, the tendency of 
to zero is characterized by the progressive domina- 
tion of its real part. The renormalized dispersion law E]^ 
(as far as the condition ^ holds) is given by the common 



(4) 



(8) 



where Ey^ = yCk+^k ^^'^ non-perturbed supercon- 
ducting dispersion law and E is specified for particular 
G^^\ j = 1, 2. Then the IRC is written down as: 

(k-ko)Vk^k»rk, 

near a nodal point ko where a nodal line crosses the Fermi 
surface. 

The low energy asymptotics of Eq|Hl correspond- 
ing to the G(2) solution, Eq|7| is: E^'' « 
(7rcl/^Po/A)ii^kln(2A/£;k), and with the related damp- 
ing r^') 



ImE(2) (^E^''^ « £;k/ln(2A/£;k), we arrive 
the condition: 



at 



Ek < Aexp 



^V^Po J 



(9) 



which defines a narrow enough vicinity of the Fermi en- 
ergy where this solution makes sense. 

Applying the same treatment to the G*-^^ solution, 
which formally defines the low energy dispersion law 



E, 



(1) 



_Ek and the damping rj^^' = ImE*^ 



(1) (4^)) 



TTcV'^pogo, we obtain the condition 
Ek > TTcV^pogQ, 



(10) 



so that this solution is valid far enough from the nodal 
points, where it provides also a correct limit of pure d- 
wave DOS. However, this solution is clearly eliminated 
near the nodal point. Thus we come to the conclusion 
that the only SCTMA solution, valid in the close vicin- 
ity of the Fermi energy, is that given by Eq[7| A physical 
consequence of vanishing DOS at e ^ for this solution 
is that the much disputed conjecture of universal elec- 
tric and thermal conductivity |25l| turns impossible. 
Nevertheless, if the validity range for the G^^-* solution, 
Eql^l is very narrow, these conductivities, as far as be- 
ing defined by the G^^^ solution, can display an apparent 
tendency to those universal values. 

Notably, the two estimates, Eas l9ll0l do not necessar- 
ily assure the overlap between the two validity regions, so 
that for ttcV^pq 3> A there can exist some intermediate 
energy range where neither of SCTMA solutions applies. 
This range roughly corresponds to the broad linewidth of 
the known impurity resonance Sres Q| where DOS cannot 
be rigorously obtained even with use of the next terms 
from GE, Eq|Sl though some plausible interpolation is 
possible between the two SCTMA asymptotics. 

Finally it is worthwhile to notice that other known 
non-perturbative solutions for d-wave disordered systems 
with DOS vanishing at e ^ as a certain power law: 
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p{e) ^ e" |j],[a|, also have to satisfy IRC since they use 
field theoretic approach, only compatible with band-like 
states. But it can be easily shown that this criterion can 
be only fulfilled for such DOS if the power is a > 1, while 
the reported values are a = 1/7 0] and a = 1 (6|. 

In fact, let the renormalized radial dispersion law (in 
the low energy limit) behave as ^ (fc — kp^ cx with 
certain v > 0, then the simplest estimate for d-wave DOS 
is 



p[e) cx £ 



cx £ 



dr] J d^S (e^ - 
drj 

(£^ — -rfY'^l'^ 



^ 1 



AvT (2 - 1,) 



that is a = 3 — 21,. In the considered field models, DOS 
defines the quasiparticle broadening Ft = u^p (^^k) i with 
a disorder parameter u. Then the criterion, Eq^ is re- 
formulated as 



» ^ P 



(a) , 



leading to the condition 3> const-^'''^^^'^^ and in the 
limit ^ ^ this is only possible if 3 — 2;^ > 1, that is 
a > 1. 

So, the above considerations essentially restrict pos- 
sible candidate solutions for quasiparticle spectrum in 
the disordered d-wave superconductor and in fact suggest 
EqUJas the only known consistent low energy solution for 
the problem. 
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